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Flow equations in generalized braneworld scenarios 
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We discuss the flow equations in the context of general braneworld cosmologies with a modified 
Friedmann equation, for either an ordinary scalar field or a Dirac-Born-Infeld tachyon as inflaton 
candidates. The 4D, Randall-Sundrum, and Gauss-Bonnet cases are compared, using the patch 
formalism which provides a unified description of these models. The inflationary dynamics is de- 
scribed by a tower of flow parameters that can be evolved in time to select a particular subset of 
points in the space of cosmological observables. We analyze the stability of the fixed points in all 
the cosmologies (our results in the 4D case already extending those in the literature). Numerical 
integration of the flow equations shows that the predictions of the Gauss-Bonnet braneworld differ 
significantly as compared to the Randall-Sundrum and 4D scenarios, whereas tachyon inflation gives 
tensor perturbations smaller than those in the presence of a normal scalar field. These results are 
extended to the realization of a noncommutative space-time preserving maximal symmetry. In this 
case the tensor-to-scalar signal is unchanged, while blue-tilted spectra are favoured. 

PACS numbers: 98.80.Cq, 04.50.+h 



I. INTRODUCTION 

The recent improvement in the determination of 
cosmological observables by the Wilkinson Microwave 
Anisotropy Probe Q, 0| and other large-scale structure 
experiments has given a boost to the search for viable 
theoretical scenarios of the early Universe. An impor- 
tant problem cosmologists still have to address is cosmic 
confusion in inflationary scenarios, whereby different un- 
derlying physics leads to the same observables; deriving 
robust conclusions from data requires an understanding 
of such model degeneracies. Viable inflation models must 
predict the quasi-invariant density perturbation spectra 
and subdominant tensor perturbations that observations 
require, but having done so it is then difficult to dis- 
criminate between different details such as the precise 
shape of the potential. Indeed, even quite radical revi- 
sions of understanding of high-energy physics, such as 
the braneworld scenario, have yet to lead to appreciably 
characteristic predictions, at least in the simplest cases 
(see Ref. Q for some reviews). 

Rather than choosing a particular class of potentials 
and calculating the preferred values in the parameter 
space, one can try to circumscribe the allowed regions 
within it through the general behaviour of the inflation- 
ary dynamics. This can be achieved, for instance, by the 
study of the consistency equations, which do not depend 
on the choice of the potential; see Ref. and references 
therein. Another possibility is to consider the evolution 
of the so-called flow parameters in terms of the number 
of e-foldings [a, la 13 • The flow equations naturally select 
a subset in the observational plane defined by the scalar 
spectral index n s and the tensor-to-scalar ratio r. Modi- 
fications to the underlying dynamics will lead to a shift in 
the location of that subset, helping to indicate whether 
a particular class of inflationary models is better able to 
produce observables in the region required by data. 



In this paper we shall explore the flow approach in sev- 
eral braneworld scenarios, where the effective Friedmann 
equation governing the cosmological brane dynamics is 
modified with respect to the 4D evolution. The simplest 
such case, the Randall-Sundrum (RS) type II model for 
a normal scalar field, was already considered in Ref. § 
and found not to introduce qualitatively different fea- 
tures. Here we shall extend the discussion to the Gauss- 
Bonnet (GB) braneworld by means of the "patch formal- 
ism" 0, which provides a unified description of a wide 
range of possible gravity models. This has the advantage 
that we can simultaneously implement general inflation- 
ary scenarios with modified Friedmann equations and ei- 
ther an ordinary or Dirac-Born-Infeld (DBI) tachyonic 
scalar field. The goal is to seek new insights into the gen- 
eral structure of braneworld and tachyon scenarios via a 
technique that might give complementary information as 
compared to model-building approaches. 

We will find that in the Gauss-Bonnet case with nor- 
mal scalar field, the flow structure in the n s -r plane is 
quite different from the 4D and RS cases because of the 
details of the flow equations. The flow structure for the 
tachyon models, here presented for the first time, gives 
independent support to other results found in literature; 
in particular, the most characteristic prediction is the 
generation of spectra with a rather low tensor-to-scalar 
ratio. The stability of the fixed points in the plane is also 
analyzed at higher order in the flow parameters. 

Finally, we shall consider how the general relativistic 
and braneworld scenarios are affected by the introduction 
of a noncommutative space-time algebra, which modifies 
the inflationary anisotropies while leaving the homoge- 
neous background (and the associated flow equations) 
untouched This skews the theoretical points in the 
n s -r plane towards the blue-tilted region, with models 
having a high tensor-to-scalar ratio suffering the largest 
shifts. 
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II. FORMALISM 



starting time). The quantity 



A. The patch universe 

We shall assume that the universe is confined into a 3- 
brane embedded in a five-dimensional bulk. Matter lives 
on the brane only, while gravitons are free to propagate 
in the bulk. This is guaranteed as long as the effective 
energy density measured by the observer is smaller than 
the mass scale of the theory, p < mf . By restricting our 
attention to the inflationary physics, we can disregard 
the contribution of the projected Weyl tensor at large 
scales. 

If the energy density on the brane is comparable to the 
potential energy adopted to stabilize the extra dimension, 
the bulk backreacts and the effective cosmological evolu- 
tion is modified by nonstandard terms. In particular, 
we can describe the primordial universe, at least in some 
finite time interval or energy patch, by the Friedmann 
equation 



H 2 



P q P q - 



(1) 



where q is constant and f3 q > is a constant factor with 
energy dimension [f3 q ] = E l ~ 2q . This equation encodes 
a number of situations, including the pure 4D radion- 
stabilized regime (q = 1, also known as the standard cos- 
mology) , the high-energy limit of the Randall- Sundrum 
braneworld (q — 2), and the high-energy limit of the 
Gauss-Bonnet scenario (q = 2/3). Defining the parame- 
ter 
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Eq. £p) can be recast as 
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= 0. 
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Deviations from the 4D case 8 = will characterize exotic 
scenarios (braneworlds or modified gravities) according 
to their magnitude and sign. 

Because of local conservation of the energy-momentum 
tensor, for a perfect fluid with pressure p the continuity 
equation is 



3H(p + p) = 0. 



(4) 



The dynamics of the patch inflationary universe has been 
extensively explored in the case of both an ordinary scalar 
inflaton <p and a DBI tachyon T; for a complete review 
of results and a list of references see Ref . |4| . It is useful 
to generically label the inflaton field as ip and introduce 
a new parameter 



for ip 
for ip 



T. 



(5a) 
(5b) 



Let Hq be the Hubble parameter evaluated at some ref- 
erence time to (later it will correspond to the integration 
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has dimension [a q ] = E^ 2 " 8 '^ 1 (= E^ 1 for ip = T and 
= E for ip = (f>), and can be absorbed in the normal- 
ization of the scalar field, so that the latter becomes 
dimensionless. In the following we use the redefinition 
ip —y ip/a q everywhere; all dropped a q factors can be 
recovered by counting the number of ip derivatives. 



B. Slow-roll parameters and flow equations 

The amount of inflation is described by the number of 
e-foldings 



N(t) = \n(at/a) 



H(t')dt' , 



(7) 



this "backward" definition measuring the number of re- 
maining e-foldings at the time t before the end of inflation 
at tf. In the flow approach, the cosmological variables 
during inflation are written as functions of N or ip; one 
can shift from one picture to the other via 



_d_ 

dN 



Hdt 



Ho 
H 



H' d 
'Wchp' 



(8) 



where a prime denotes differentiation with respect to ip. 
We define the flow parameters as 



A 



_ d\nH 

= dN 

H 0/ 



H- 



Ho 
H 



EL 



Ho 
H 



HV +1 \H'] 



£-1 



t > 1 



(9a) 
(9b) 



where (n) indicates the n-th ip derivative. Usually in 
literature they are dubbed Ajj. Note that this definition 
can be obtained from the potential slow-roll (SR) tower 
of Ref. Q with the substitutions 4> => ip, V ^ H, and 
q — > 1 + 8. In four dimensions and with an ordinary scalar 
field this definition coincides with that of Ref. |(| . 

It is also convenient to consider the following quantity: 



a = 2Ai - (4 + 8 + 6)e. (10) 
The set of evolution equations reads, from Eqs. (jSJl, 
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and (TUJl, 



de 

In 

da 

In 

dN 



e[(r + (2 + 6)e] , (11a) 

2A 2 -(5 + 29 + 0)ea 

-(A + + 0)(2, + + 0)e 2 , (lib) 

\ e+1 +l\ e {(l-l)a 

+ [2£0-(4 + + 0)}e}, £>2, (11c) 

where 20 = 2 + 0-6. 

As shown in Ref. [7J , the flow parameters and equations 
do not encode the inflationary dynamics, since they are 
only a set of identities for the Hubble rate H as a func- 
tion of N. In fact, by definition the next-to-lowest-order 
SR terms appear with higher powers, and one can ap- 
proximate the dynamics through a power truncation of 
the traditional SR tower. In the case of the flow param- 
eters, Eq. this is achieved by imposing a constraint 
such as = and H^> ^ 0, for some maximum 

n. However, the method still provides an algorithm for 
generating inflationary models in the space of the cosmo- 
logical observables, corresponding to a Taylor expansion 

of#(VO0- 

For reference, the flow parameters are related to the 
Hubble-slow-roll parameters, defined by |4l.lT(i| 
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1/r, 



n > 1, 



via the relation 



(») 



(12) 



(13) 



For a comparison between the two towers from an obser- 
vational point of view, see Ref. . 



C. Attractors and observables 

It is easy to find the attractors of the flow and project 
them on the a-e plane. One is the line corresponding to 
de Sitter (dS) fixed points 

e* = 0, (14a) 
a* = const, (14b) 
A; = 0, £>2, (14c) 

and the other is the power-law inflationary attractor 

e* = const , (15a) 
a* = -(2 + 0)e, (15b) 



a; 



(l + 0)(2 + 0)e 2 , 
2 + (1 + £)§} X* e e* 



£>2. 



(15c) 
(15d) 



Note that the slope of the line Eq. (|15bjl depends on the 
braneworld model one is considering, but not on the type 
of scalar field on the brane. 

To first-order in slow roll the cosmological observables 
are 



e 



1 



(16) 



where r = A 2 /A 2 is the ratio between the tensor and 
scalar perturbation amplitudes and n s is the scalar spec- 
tral index. The coefficient ( q is equal to 1 in the 4D 
and GB cases while C2 = 2/3 in the RS scenario 
|13| . The expression for r is valid in all cases considered, 
while the expression for n s is for the commutative geom- 
etry case only and will be modified later to address the 
noncommutative case. 

The stability analysis of the above fixed points against 
linear perturbations in the flow parameters involves an 
infinite-dimensional parameter space. However, one can 
truncate the flow tower at ^-th order and study the sys- 
tem order-by-order, seeking convergence of dynamical 
properties. 

Define the vector X = (e, a, X 2 , . . . , Xe) 1 . The pertur- 
bation equation reads 



dSX 

In 



MSX, 



(17) 



where SX = (de, Sa, . . . )* encodes the linear perturba- 
tions and M is a (£ + 1) x (£ + 1) matrix whose elements 
m ij (*! 3 = 0,...£) are evaluated at a fixed point. For 
i < 2 (the first two rows) one has 

m 00 - 2(2 + 6>)e*, m i = e*, 
m 10 = -[(5 + 20 + 0)<r* +2(4 + + 0)(3 + t 



mix = -(5 + 20 + 0)t*, mi2=2, 
m*j =0, j > i + 2, 

while for i > 2 



m i0 
ma 
ma 

mij 



-[2i9-(A + e + 9)]X*, 
l{(i-l)a* + [2i6-(4 



(18a) 

(18b) 
(18c) 
(18d) 



(18e) 
(18f) 



(18g) 

1, (18h) 
0, 1 < j < i, j>i + 2. (18i) 



The stability condition requires that the real part of 
all the eigenvalues 7 of the matrix M is nonnegative. 1 



In the following the operation 7 — » SR(7) is understood. 
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The sign for the stability depends on the convention set 
by Eq. (J7J; nonpositive eigenvalues correspond to fixed 
points that are stable in the past. 

For the dS fixed points, the Jordan equation 
det(M — 7) = is simply 



7 



' 7 



(£-1)<j* 



7 



0. 



(19) 



and stability is guaranteed for a* > (blue-tilted spec- 
tra). 

The power-law case is more complicated since one can- 
not easily diagonalize the matrix. The only notable ex- 
ception is the Gauss-Bonnet braneworld with a normal 
scalar field. In that case, A| — for £ > 2 and the 
eigenvalue equation reads 



r 



201 ¥ -7 



(£ - l)e* 



7 



= 0. (20) 



The first two eigenvalues are equal in absolute value and 
opposite in sign. Therefore the power-law fixed points 
in GB (^-inflation are a repcllor in both backward and 
forward integration. 

In other cases of interest a full analysis is required; 
moreover, one should check the stability in both back- 
ward and forward time integration and along a "suffi- 
cient" number of directions in the parameter space. If 
all the eigenvalues have the same sign, then there is an 
attractor either in the past or in the future. Conversely, 
if some of them have relative opposite sign, there is no 
attractor. The truncation order is also important, since 
higher-order terms can change the sign of the eigenvalues. 

To determine the minimum level at which the stability 
analysis is reliable, we have solved the Jordan equation at 
the power-law fixed points from first order (2x2 matrix) 
to 10 th order in the flow parameters (£ = 10), and checked 
that the real parts of the eigenvalues have different signs 
at higher orders. Figure ^ shows the eigenvalues of the 
4D and RS models with a normal scalar field and the 
GB tachyon case. The GB (j> case is as described above, 
while the other tachyon plots are quite similar to that 
presented here. 

Qualitatively these results do not depend on the actual 
value of e*; the only assumption is e* > 0. As regards 
the level of approximation, for the normal scalar field 
and the 4D tachyon the first order is already sufficient to 
give the result holding at all higher orders, while in the 
RS and GB tachyon cases all eigenvalues are negative up 
to 2 nd and 3 rd order, respectively. At higher orders one 
eigenvalue becomes and remains positive. Note that the 
degeneracy of one of the eigenvalues (the lowest one in 
the RS (j> and GB T cases) is lifted every two orders (at 
even orders in the cited example). 

The integration of the flow equations, performed at 5 th 
order (see below), confirms these results. Therefore one 
can conclude that the power-law fixed points are unstable 
in both time directions, for all the present cosmological 
models, and at any (high) order in the flow truncation. 
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FIG. 1: Real part of the eigenvalues of the perturbation 
matrix for power-law fixed points, from 1 st to 10 th order in 
the flow parameters. From top to bottom, the 4D <f>, RS 0, 
and GB T models are shown for e* = 0.1. 



The latter has been shown to be consistent in all cases 
from 4 th order on. 



III. NUMERICAL INTEGRATION OF THE 
FLOW EQUATIONS 

In order to study the evolution of the flow parameters, 
we have to establish a truncation level for the flow hi- 
erarchy and set suitable ranges of initial conditions for 
the parameters themselves within which they are ran- 
domly selected. Early flow papers solved the flow hier- 
archy numerically but a more efficient approach [T^ | 
exploits an analytic solution to the flow equations found 



5 



in Ref. 0: 



H(iP) = H Q 



M+l 

1=1 



(21) 



where Hq = H(to) is evaluated at the initial integration 
time to; without loss of generality we choose ip{to) = 0. 
This solution continues to be valid in all cases considered 
here. With the above normalization for the scalar field 
the coefficients Ac are dimensionless: 



A.= 



(-1) £ A, 



e>i, 



(22) 



V>=o 



where A^_i = o = 1- The choice of the — sign, coming from 
the (H'Y -1 term in Eq. I|9bj) . is a convention determined 
by the rolling direction of the scalar field, which in this 
case is such that ip > 0: 



H 



0/2 



(23) 



Note the coefficients Eq. (|22|l have a different normal- 
ization with respect to Ref. [l4|. due to dimensionless 
factors such as ^/4n absorbed into the definition of if). 

To obtain the full evolution of the flow parameters, one 
needs only numerically integrate the relation between N 
and if), which from Eq. (jSJ) is 



dN 




(24) 



For each model we used the same sequence of 40,000 ini- 
tial conditions chosen within the ranges [|J 



eo 

(TO 
A2,0 

A3.0 



[0,0.8], 
[-0.5,0.5] , 
[-0.05,0.05] , 
[-0.005,0.005] 



(25) 



A, 



6.0 



= 0. 



For an implementation at 5 th order in slow-roll the value 
of the last parameter truncates the series and closes the 
hierarchy. Initial conditions which do not generate either 
a sufficient amount of inflation or a well-defined asymp- 
totic behaviour are rejected. Conversely, when enough 
inflation is realized, then either 

1. The evolution proceeds until the end of inflation, 
e(U) = 1. Then the flow equations are integrated 
backwards for a suitable number of e- foldings where 
the observables n s and r are evaluated and plotted, 
or 

2. A late-time attractor is reached and the observables 
are read off at that point. 



TABLE I: Population counts of forward integrated, backward 
integrated, and attractor points in the plots of Fig. as 
percentages of a set of 40,000 entries. 



% 


No > 50 


iV < 50 


Attractor 




4D 


0.2 


6.8 


93.0 


<A 


RS 


0.1 


4.5 


95.4 




GB 


0.5 


16.7 


82.8 




4D 


0.1 


4.6 


95.3 


T 


RS 


0.1 


3.9 


96.0 




GB 


0.1 


5.7 


94.2 



The reference number of e-folds for the back- 
ward/forward integration can be cither fixed or randomly 
chosen within a range. Here we take No = 50, although 
the results will not depend significantly upon its precise 
value 0. Also, we can distinguish points evolving to 
£f = 1 in two subcases: 

la. More than 50 e-foldings of inflation are obtained 
from the initial point, meaning that the location 
where the observables are read off was reached by 
forwards integration from the initial condition. 

lb. Less than 50 e-foldings are obtained from the ini- 
tial point, so that the point corresponding to the 
observables is obtained by integrating backwards in 
time from the initial condition. 

The following results were all obtained via the ap- 
proach described above. Additionally we checked that 
integrating the full flow hierarchy gives the same results 
point-by-point, which is a useful cross-check of the nu- 
merical implementation. 



IV. RESULTS 

We divide our results into the (standard) commutative 
case and the noncommutative case. In the latter the flow 
equation evolution is unchanged, but the expression for 
the spectral index Eq. Ijltif) is altered. 



A. The commutative case 

Figure |21 shows the flow points in the n s -r plane for 
an ordinary and a tachyonic scalar field in the 4D, RS, 
and GB scenarios. The classification of points is model 
dependent and summarized in Table [I] In all cases the 
attractor points are the dominant class, and hence the 
blue-tilted region (n s > 1) is more populated than the 
red-tilted one. 

The points that finish inflation with 50 e-foldings or 
more integrating forwards in time are a very small group 
of the three considered categories. Most of the points 
that finish inflation do so with less than 50 e-foldings. 



6 




.3 i 



(c) GB, 



FIG. 2: Distribution of commutative inflationary observables in the n B -r plane. Left column: plots for an ordinary scalar field 
4> in the a) 4D, b) RS, and c) GB cases. Right column: plots for a DBI scalar field T in the d) 4D, e) RS, and f) GB cases. 
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Note that what one does with the integration backwards 
is to find the point in parameter space where the initial 
condition should have been in order to give 50 e-foldings 
or more. This does not necessarily correspond to the 
process that gave rise to the physical initial conditions 
for inflation to start. 

Figure|2roughly shows that the swathe near the power- 
law attractor decreases for increasing 9, so that its ver- 
tical extension decreases from GB to 4D to RS. This is 
clear from Eq. (|lla|l . since the second positive term in 
square brackets becomes more and more important. This 
seems to indicate faster flow, which would naturally give 
less height to the swathe as the natural flow is towards 
the r = line. Note that the swathe corresponds to the 
flow being integrated backwards in time from the initial 
point in order to generate 50 e-foldings. 

Comparing the percentages, it is clear that in the GB 
(j> case there is an increasing number of backward inte- 
grated points, eating off about 10% of the population 
of attractor points. Despite this fact, the swathe loses 
coherence and points spread in a large region in the ob- 
servable plane, making the above descriptions rather ap- 
proximate. Although there is no particular feature in 
the GB power-law attractor, the plot points do not lie 
close to it as in the other cases. This feature is a direct 
consequence of the flow equations. Equation shows 
that the relevant parameter governing the cosmological 
evolution is 9 rather than q; then the GB flow (6 = —1) 
would act in the opposite direction with respect to the RS 
(9 = 1) case, the 4D one being intermediate between the 
two. More precisely, while in RS the points are more con- 
centrated at the region around r — 0, n s = 1, it is natural 
to expect more a widespread point distribution for GB. 
This does not necessarily imply that the GB braneworld 
is more severely constrained than the RS one. 

In the GB system, contrary to what happens in the 
others, all equations can be written in terms of A2 and 
integrals of it. Since 



A 2 = 



dN 



+ (l + 0)eAx, 



for 9 = 9 = — 1 = — 9 one has 
A x = J X 2 dN, 



(26) 



(27) 



A, = I + (£ - l)\i\t]dN, £>2. (28) 

This property makes the numerical integration of the GB 
equation peculiar relative to the other models. Also, as 
shown in Sec. Ill Cl the extra directions orthogonal to the 
n s — r plane completely decouple in the dS and power-law 
cases (Xg ~ for I > 2). Therefore the two-dimensional 
analysis presented above is expected to encode all the 
relevant dynamical features. 

In the tachyon case, the percentages do not change 
appreciably in different braneworlds, nor with respect to 
the ordinary scalar situation. This is because Eq. (JSJ) is 



insensitive to the type of braneworld in the tachyon case 
and the information on the extra dimension is carried 
only by the parameter Ai. Therefore the tachyon plots 
are very similar to one another. 

Also, in this case plot points are squeezed towards the 
dS attractor and prefer lower values of r. In order to 
understand why, let us define q\ = a 2 + r 2 ; the difference 
of the (squared) radii for a fixed braneworld {Cg,#} is, 
from Eq. l(TU|) . 



A 2 = 



q\ = (2-6>)[4Ai - (lO + 30)e]e. 



(29) 



For a given input set {e, Ai, • • • } of parameters, A 2 > 
in any case of interest considered so far and for the main 
part of initial conditions, which explains why the tachyon 
points are concentrated near the scale-invariant origin 
of the a-r plane, while <f> points are more widespread 
through the plane. We verified numerically that A 2 < 
in less than 10% of the cases. A similar argument, with 
the SR parameters now depending on the scalar field, can 
be applied for the large-field model plots of Ref. [l5| . 

Note that these results are completely determined by 
the choice of the patch parameter 9. The parameter £ g 
would need to be extracted from the two-point correla- 
tion function of the graviton zero-mode within the orig- 
inal gravitational theory (3-brane in Einstein or Gauss- 
Bonnet gravity). However, one can redraw the plots in 
terms of the quantity £ 9 r = 0(r), which differs from r 
only in the RS scenario. In this case, the only effect 
is a further squeezing (by a factor 2/3) of those points 
towards the dS attractor. 



B. The noncommutative case 

Beside extra dimensions, one can take into account a 
modification of the space-time geometry at the quantum 
level. In particular, a space-time uncertainty principle is 
believed to be a universal property of high-energy theo- 
ries motivated by strings [l6| . This can be realized by 
a noncommutative algebra which preserves homogeneity 
and isotropy in the cosmological context [jj : 



[r, x] = il 2 s 



(30) 



where r = J a dt is a redefinition of time, x is a spatial 
coordinate on the brane, and l s is the fundamental string 
length. In the simplest version of this model, moduli 
fields are fixed so that the extra dimension commutes 
with all the other directions. 

The presence of a particular scale further breaks scale 
invariance and the perturbation spectra are affected ac- 
cordingly. The detailed consequences of this ansatz were 
explored elsewhere @,llMll3- Here we just recall that at 
large scales (corresponding to the strongly noncommu- 
tative limit) there are basically two different implemen- 
tations for noncommutativity ("class 1" and "class 2"), 
depending on whether the Friedmann-Robertson- Walker 
2-sphere is factored out from the measure of the action 
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for the cosmological perturbations or not. In both cases 
the background, flow equations, percentage populations 
of Table [I] and stability analysis are not modified, and 
the new ingredient is encoded only in the expression for 
the scalar spectral index: 

n s — 1 w a + ce, (31) 

where a is defined in Eq. (|l(Jfl and the constant c is c = 
in the commutative case, c = 2 in the class 2 noncommu- 
tative model, and c = 6 in the class 1 model. Since the 
scalar and tensor amplitudes are multiplied by a same 
extra factor depending on l s , their ratio r is unchanged. 
Therefore the points in the upper region of the n s -r plane 
are shifted towards the right (blue-tilted spectra), while 
those lying on the r = line are unaffected. As one can 
see in Fig. |3 for c = 6 the swathes point sharply right- 
wards instead of leftwards as in the commutative case. 
For c = 2 the swathes are almost vertical. Thus the ef- 
fect of noncommutativity is a skewing in the plane whose 
magnitude depends on the value of the parameter c |15| . 

Note that since the noncommutative correction to the 
spectral index is positive, the stability condition for the 
dS fixed points is again satisfied for blue-tilted spectra. 

For brevity we have not shown tachyon plots in the 
noncommutative case, but their shape can readily be pic- 
tured by applying the same skewing to the points shown 
in Fig. 

C. Relation to observations 

We conclude by comparing the above flow plots with 
observations, with an analysis similar to that of Peiris 
et al. 0] and Makarov ^l| for the standard cosmology. 
Present observations are highly restrictive in the n s -r 
plane, with only a small region near the origin still al- 
lowed in the plots we showed (very roughly, the allowed 
region corresponds to r < 0.03 and 0.95 < n s < 1.1, but 
one must allow for those parameters being correlated). 
The regions where our plots show major differences are 
thus already excluded. We take the 95% likelihood bound 
of Ref. , conservatively choosing the wider region ob- 
tained without including Lya-forest data, and compute 
in each case the fraction of points lying within that re- 
gion (which we approximate by a suitably-oriented el- 
lipse, noting that their r is defined to be 16 times ours). 
We do not impose a constraint from the spectral index 
running. 

The percentage of points inside the allowed region is 
shown in Table ITT1 In all cases the percentages are quite 
small, indicating that observations have already chopped 
off a substantial part of inflationary model space. Bear- 
ing in mind that the majority of the points are late-time 
attractor points with n s > 1 and negligible r, the strong 
constraint n s < 1.04 at r = pl a Y s a major role in 
ensuring the percentages are so small. If this limit were 
further strengthened in the future to exclude blue-tilted 
models, all those points would be lost. 



TABLE II: The percentage of points lying within the 
currently-allowed region for each of our cosmologies. 
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0.6 
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There are a couple of trends evident in the numbers. 
Firstly, the tachyon cases tend to give higher percentages, 
due to their suppression of the tensor ratio r. Secondly, 
for the normal scalar field, increasing the noncommuta- 
tivity parameter c increases the percentages as it skews 
the swathes across to where the observations lie. Nev- 
ertheless, every cosmology proves capable of generating 
models within the allowed region, which prevents any 
useful conclusions from being drawn. 

V. CONCLUSIONS 

In this paper we have analyzed braneworld and tachyon 
scenarios in terms of the flow evolution equations. The 
equations of motion and the stability of the associated 
fixed points have been considered in detail in the general 
relativistic, Randall-Sundrum, and Gauss Bonnet cases 
at high order in the flow parameters. Each model gen- 
erates different predictions for the main cosmological ob- 
servables (scalar spectral index n s and tensor-to-scalar 
ratio r). In particular, Gauss-Bonnet gravity deeply 
modifies the flow structure in the parameter space for an 
ordinary inflaton field. In the cosmological tachyon case, 
the theoretical points in the n s -r plane are pressed to- 
wards scale invariance but the characteristic flow swathe 
does not lose coherence. We also discussed the imprint of 
a noncommutative geometry on the inflationary observ- 
ables. 

Some trends remain to be explored. In the braneworld 
inflationary scenarios the Weyl tensor is negligible at 
large scales and we have consistently neglected its con- 
tribution to the observables r and n s . However, at small 
scales it can play an important role, as well as in the 
determination of other observables such as the bispcc- 
trum even at long wavelengths. Although we expect the 
flow approach to be unmodified by bulk physics, further 
investigation might clarify this point. 

The flow equations approach is only one possible way 
in which one can generate an ensemble of inflationary 
models, corresponding to a Taylor expansion of H(ijj). 
In Ref. |l4j , other methods were implemented in the gen- 
eral relativistic case to study the robustness of the flow 
equations structure, showing that there are significant 
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FIG. 3: Distribution of noncommutative inflationary observables in the n s -r plane for an ordinary scalar field (j>. Left column: 
plots for the class 1 model c = 6 in the a) 4D, b) RS, and c) GB cases. Right column: plots for the class 2 model c = 2 in the 
d) 4D, e) RS, and f) GB cases. 
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variations in outcome between methods. It may be in- 
teresting to explore these differences in the braneworld 
and tachyon cases, too. 
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